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We investigate the influence of domain walls on the vortex dynamics in superconductors 
with multi-component order parameters. We show that, due to their complex structure do- 



r-{ • main walls can carry vortices with fractional flux quanta. The decay of conventional vortices 

j/3 ' into fractional ones on domain walls is examined. This decay presents an extraordinarily 

strong pinning mechanism for vortices and turns domain walls occupied with pinned frac- 
C^ , tional vortices into efflcient barriers for the vortex motion. Therefore, domain walls can 

act as fences for the flux flow, preventing the decay of the remnant magnetic flux enclosed 
by them. Furthermore, the consequences of this property of domain walls on the vortex 
i-rt ', dynamics are discussed in connection with observed noise in the hysteresis cycle, using the 

r-j ■ Bean model of the critical vortex state. Based on this picture experimental data in the 

Q ' unconventional superconductors UPta, Ui_a:Tha;Bei3 and Sr2Ru04 are interpreted 

o : 

§1. Introduction 

> 

^ ■ The penetration of magnetic fields into type-II superconductors as flux lines 

^^ . (vortices) yields many complex phenomena. In recent years the physics of vortices 

(<— ^ I has become an imnprtant subject, particularly in connection with high-temperature 

T-H ■ superconductivity. EJ) The physics of "vortex matter" is crucial for a large number 

On , of applications of superconductivity involving high currents and fields. Vortices 

-■.^ I limit the technological capability of superconductors, since their motion generates 

C^ ' dissipation. Therefore one of the central issues is the pinning of vortices at defects 

such as impurities, lattice dislocations and twin boundaries. The critical current, 
, ^ the limit of dissipation-free transport, depends on the character and strength of the 

f^ ■ pinning potentials created by these defects. The effect of pinning was shown to 

O ■ crucially depend on the variijus phases of vortex matter, and the resulting problems 

are naturally very complex. EP 

In this paper we would like to discuss a pinning phenomenon in unconventional 
K^ , superconductors which has a character different from the usual pinning at crystal 

J^ I defects. This study is motivated by recent experiments on heavy fermion comnounds 

UPta and Ui-^Th^^Beis as well as on the transition metal oxide Sr2Ru04B''l5''Q''EP. 
All these superconductors have a comparatively low transition temperature of the 
order of IK. Fluctuation effects do not play an important role in these systems. 
However, interesting vortex physics is introduced because the superconducting order 
parameter of these systems has more than one component. This yields more degrees 
of freedom in forming topologically stable defects of the order parameter, a well- 
typeset using VTVT^.sty <ver.0.8> 
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known feature in superfluid ^He physics. We would first like to review some of the 
basic experimental facts concerning the mixed state of these three superconductors, 
which will be relevant for our theoretical discussion. 

One aspect of the mixed state related to vortex pinning is the slow motion of 
vortices close to the critical state, known as flux creep. This phenomenon is observed 
as the slow decrease of the remnant magnetization after the superconductor has been 
exposed for some time to a magnetic field considerably larger than the lower critical 
field, Hci- The remnant magnetization that exists after turning off the external 
field originates from the pinning of vortices in the material. The motion of these 
vortices out of the sample, generally governed by thermally activated crossing of 
pinning potential barriers, leads to the slow decay of the magnetization. HP (Note that 
quantum phenomena, macroscopic quantum tunneling, can also play an important 
role at sufficiently low temperature. EP'S-') The decay (creep) rate is determined by 
the temperature, the pinning properties, and the vortex rnatter state as described by 
Kim and AndersonlZP'EP (see also Geshkenbein and LarkintP). A simple and intuitive 
theory of the critical vortex state in a material with strong pinning effects was given 
by Bean. [3 It describes the characteristic profile of the vortex distribution depending 
on the history of the external applied fields and the critical current. eP 

In the two heavy Fermion superconductors UPta and Uo.9725Tho.o275Bei3 (0.017 < 
X < 0.45), an anomalous temperature dependence of the creep rate was observed by 
Mota's group. |)'& UPts and Ui_:i.Th2;Bei3 (0.017 < x < 0.45) both exhibit two 
consecutive superconducting phase transitions. Flux relaxation measurements show 
a rapid drop of the creep rate, essentially to zero, immediately below the second 
transition in both system. lP'H^ A similar sudden drop in the creep rate was reported 
for Sr2Ru04 at the rather low temperature T* ~ 50 mK.o' In contrast to the former 
two examples, no sign of an additional transition at this temperature is observed in 
any other properties to this time. The transition in the creep rate indicates the onset 
of a new efficient pinning mechanism which inhibits the motion of vortices from the 
interior to the surface of the sample. However, it has been noticed in experiments 
that after a long waiting period some creep recovers. B''q' 

This effect is readily explained if we assume that fence-like structures exist which 
prevent the passage of vortices (so that the vortices cannot leave the sample) . Ex- 
periments indicate that these fences are activated below a certain transition temper- 
ature within the superconducting phase. The effect of the fences is rather different 
to standard pinning. In particular, their influence on the critical current is weak, as 
discussed below. This is consistent with hysteresis measurements on UPts in a slowly 
oscillating magnetic field. B> While the flux creep disappears in the low-temperature 
phase, the critical current is continuously increasing and shows no anomaly upon 
entering the low temperature phase. B^'H' Furthermore, below the second transition 
Rosenbaum's group found that the magnetization curves of the hysteresis cycle dis- 
play strong variations from cycle to cycle in the region close to Hd , while in the high 
field part of the cycle (connected with the critical current of the superconductor) no 
unusual behavior is seenEP. This "noise" can be understood as being induced by the 
weakly mobile fences inhibiting the free flow of vortices. In every cycle the position 
of the fences changes slightly, and the magnetization process for small vortex con- 
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centrations, i.e. for fields close to Hd, is modified. 13^ We will show that a simulation 
with a simple model incorporating this effect gives rise to the features observed in 
these experiments. 

What is the origin of these fence-like structures? In unconventional supercon- 
ductors such structures occur if the superconducting phase is degenerate, e.g., for 
time-reversal symmetry breaking states which have at least a two-fold degeneracy. 
Degenerate states can appear as domains in the superconductor accompanied by 
domain walls as topological structures of the order paranietenEiP Various properties 
of such domain walls have been studied by several groups. HIP 'Ej)'E3* Domain walls 
in time-reversal symmetry breaking superconductors npasess interesting magnetic 
properties and carry chiral quasiparticle bound states. O'H'Ej) The question arises 
whether or not they serve as barriers to the vortex motion. The domain wall is 
naturally accompanied by a slight local suppression of the order parameter which 
tends to trap vortices. In this case the pinning strength is rather weak and should 
not have much influence on the flux motion. However, strong pinning can arise from 
properties of the internal structure of domain walls. Under certain circumstances the 
domain wall has, in addition to a stable structure, metastable structures, or even two 
degenerate stable states. Hll' These states may form domains on the domain wall that 
are separated by line defects. It has been shown that such lines carry a magnetic fhix 
which is an arbitrary fraction of a standard flux quantum in a superconductor. 113 
We will show that an ordi nar y vortex placed on such a domain wall can decay into 
two fractional vortices. Oj)' 113 Since these fractional vortices can only exist on the 
domain wall, they are strongly pinned. They repel other approaching vortices, and 
in this way the domain wall indeed acts as a strong barrier. A high density of vortices 
close to the domain wall can destroy the state which carries fractional vortices, and 
the pinning effect disappears. This is equivalent to having a barrier "height" for the 
vortex pinning. It is important to note that this pinning property of the domain wall 
does not alter the bulk pinning due to impurities within each domain. This type of 
impurity-based pinning is responsible for determining the critical current. However, 
the domain walls play an important role in the overall flux motion. They are the 
origin of an intrinsic pinning location created by the superconducting state itself, in 
contrast to the extrinsic pinning due to material defects. Note, however, that the 
domain walls themselves are pinned at lattice defects reducing their mobility. 

In this article we first discuss properties of domain walls in the most simple case 
of a time-reversal symmetry breaking superconducting state using the Ginzburg- 
Landau theory. Then we use the domain walls as barriers in a model based on the 
Bean's theory and consider the modification of flux motion properties due to these 
barriers. 

§2. Properties of a domain wall 

In this section we investigate the structure of a domain wall in a time-reversal 
symmetry breaking superconducting state as may be realized in the low-temperature 
phases of UPts and Ui-j^ThajBeis, and at the onset of superconductivity in Sr2Ru04. 
We do not include the aspects of the double transition related to the former two, as 
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it complicates the discussion considerably without leading to further insight. 

2.1. Ginzburg- Landau theory 

For the following discussion we introduce the simplest possible model that con- 
tains all the relevant features in order to discuss the physics of domain walls in an 
unconventional superconductor. We consider a system with tetragonal crystal sym- 
metry based on the point group D^fi- The superconducting order parameter shall 
belong to the two-dimensional representation Eg (even parity) or Eu (odd parity), 
yielding the following expansion for the gap functions (2x2 gap matrix with the 
notation A = iayip for even and A = i{d ■ cr)(Ty for odd parity): 

1p{k) = {r]xVzVa: +riyVzVy)/{vlvl), 

(2-1) 
d{k) = 7.{rixVx + VyVy)/{vl)- 

Here T} = {rix,i]y) represent the two-dimensional complex order parameter, z denotes 
the unit vector along z-direction, Vi denotes the components of the Fermi velocity, 
and (...) is the average over the Fermi surface. As mentioned above, we concentrate 
on the time-reversal symmetry breaking states of the form r/ oc (1, ±i) yielding the 
gap functions 



V'(fc) = VoVzivx ± ivy)/{vlv'^,) or d{k) = z{vx ± ivy)/{vl). (2-2) 

It is convenient to transform the order parameter into the basis set 'q± = {rjx =F 
ir]y)/y/2. For simplicity we also use the phenomenological parameters as determined 
by weak coupling theory and assume that v^ does not depend upon kx or ky. The 
general GL free energy then has the form 



F= fd^ [a(|r/+|2 + |r/_|2)+6{(|r/+|4 + |r/_|4)+4|??+|2|r/_ 



+ i((D_r/+)*(D+r?_) + z.(i^+r/+)*(Z)_7?_) + c.c.)} + ^(V x Af]. 

(2-3) 
Here a = ao{T — Tc) and b and k are the standard coefficients derived from the 
weak coupling theory. In the gradient terms, we use the gauge invariant derivatives 
D = 'V—i^A, with D± = DxUDy and 7 = 2e/hc = 2tt/(Pq (^0 is the flux quantum). 
The parameter u denotes the deviation of the Fermi surface or its density of states 
from cylindrical symmetry around the z-axis due to the tetragonal crystal field 

{vt) - 3{vM) , ^ 
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where z^ = for a cylindrical symmetric Fermi surface. In this form the GL the- 
ory describes two degenerate superconducting states, (?7+,?7_) = riQ{T){l,Q) and 
r/o(r)(0,l) with 

ril{T) = '-^^^. (2-5) 

for T < Tc. 

In the following we analyze the properties of domain walls between the two 
degenerate states. We assume that they are infinite and planar such that we can 
separate the coordinates into components parallel and perpendicular to the normal 
vector n of the domain wall. For simplicity we ignore the z-direction, assuming that 
n _L z. Therefore it is convenient to rewrite the free energy in these coordinates: 
(x,y) — > {x',y') = (r|j,r_|_) and D± = e^^^D'^. Simultaneously we transform the 
order parameter ri± = e rj'^ so that the gradient terms become 

and in the homogeneous part only the third term among the fourth order terms is 
modified to 

6z.e^^V-!N+ + c.c. (2-7) 

From this point we omit the primes for coordinates and order parameters and take 
9 as the angle of the normal vector n relative to the crystal x-axis. Thus, the 
a;(y)-coordinate is now always parallel (perpendicular) to the normal vector. In this 
representation the 9 dependence clearly does not appear, if the system is cylindrically 
symmetric, i.e., u = 0. 

2.2. Structure of the domain wall 

We now turn to the domain wall and consider the situation in which for x — > ±00 
the superconducting states r]± is realized and around x = a smooth change between 
the two states occurs within a finite length scale (Fig. ||). This situation is similar 




X 

Fig. 1. Schematic behavior of the order parameter components at the domain waU. 

to a Josephson junction between two superconductors with the state 77+ and 7/_, 
respectively. Indeed the phase difference between the two states plays a similarly 
important role for the physical properties of the domain wall. The structure of the 
domain wall is obtained by a variational minimization of the free energy J^. For 
the sake of transparency we simplify the treatment by introducing the following 
approximate variational ansatz for the order parameter: 
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7?+ = r/oe*'^+ cos X and ??_ = ryoe*'^" sinx- (2-8) 

This leads to the boundary conditions 

{0 X ^ +00, 

(2.9) 
- x^ -oo 

We assume spatial dependence only for x and leave the phase difference between the 
two sides a = (/>+ — 0_ constant. The effective free energy becomes 

:F = Jd^r [a7]^ + 6r/^ + ^(1 + 1/ cos(2a + 46*)) sin^ 2x 

+K77o{|-DcosxP + |I?sinxP + [e"*"(L'_ cosx)*(-D+ sinx) 

+i/e-^("+^^)(L'+cosx)*(-D- sinx) + c.c]/2} + — (V x Af]. 

Svr 

(2-10) 
We neglect contributions from Ay and vary the free energy with respect to X; 

a^x = |sin4x-45+7.4,.(9,x) + 2C7+(sin2x92x + cos2x[(a,x)' + 7'4]), (2-11) 



and Ax, 



with 



^x = ,,^t. ^ y 2-12 

7(l + C+sm2x) 



Q =j^{l + i^ cos{2a + 4:9)}, (2-13) 



1 

eg 
C± = [cos a±i' cos{a + 46*)] /2, 

S± = [sin a iu sin(a + 26)] /2, 



and ^0 = '^f^^'Ho- The third term on the right-hand side of Eq. ( 2-ll| ) and the sin 2% 



term in the denominator of Ax complicate the solution. We find, however, that the 
influence of these terms is small at larger distances from the domain wall, as one 
can see in approximating x{^) ~ e~^^ in the limit x -^ +cxd. With this ansatz, Eq. 
( |2-11| ) becomes 

(1 - Sl)K^e-'^'' = ge"^^ + 0(e-2^^). (2-14) 

Thus the second term decays much faster than the first one. Therefore, the width of 
the domain wall is largely determined by the first term, and the second term enters 
only into the domain wall energy. 
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Fig. 2. The domain wall energies as functions of the phase difference a = 4>+ 
9 — 0, 7r/8 and 7r/4. The anisotropy parameter is chosen f — —0.5. 



for angles 



We solve Eq. (2-11) neglecting the third term and the sin 2% term in the denom- 
inator of Ax, so that it takes the form of an ordinary Sine-Gordon equation, 



^2 Q 

9xX = -rsmAx 



(2-15) 



where Q = Q/{1 — S'^) for which we easily find the kink solution 



xix) = arctan(e 



-x^/Q\ 



(2-16) 



This is used to determine A^ through Eq. ( |2-12| ) and both Ax and x are inserted 
into the free energy and integrated, so that we obtain for the variational domain wall 
energy 



f{a,e) = ^r]lJQ{l 



C:7r Si 



C. 



1 



l-C! 



=arctan 



1 + C+ 



(2-17) 



as a function of a and 9. 

From this point we fix z^ < 0. For sufficiently large \iy\ the domain wall energy 
develops two local minima as a function of the relative phase a. This is illustrated 
in Fig. 2 for angles 9 = 0, tt/S and 7r/4. For the most stable orientation 9 = these 
minima are at a = (stable state) and a = ir (metastable state). For 9 -^ 7r/4 we 
find that the stable and metastable state approach in energy and become degenerate 
at exactly 9 = 7r/4. The presence of stable and metastable states is important when 
we discuss vortex states on the domain wall. 

§3. Vortices on the domain wall 



We now investigate local modifications of the domain wall structure which give 
rise to vortices. This problem has been previously considered for the case of the 
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degenerate domain wall state. Ej-' Here we extend the discussion to the general situ- 
ation. 

3.1. Fractional flux lines 

As mentioned in the Introduction, domain walls act as pinning regions for the 
vortices (core pinning) due to the locally diminished condensation energy. The re- 
sulting pinning effect is, however, rather small, since the domain wall only provides a 
shallow attractive potential. In the following we would like to show that a change of 
the vortex structure can lead to considerable strengthening of the pinning potential. 

To some extent the domain wall can be viewed as a planar weak link between 
two superconductors with order parameters slightly interpenetrating, as described 



by Eq. (2-8). We use this form for the order parameter to express the current which 
is obtained from the variation of the free energy with respect to the vector potential, 



jx = — [2 cos^ X'«'a;+ + 2 siu^ x^o;- - sm2x{iux+ +Ux-)C+ + {uy+ + Uy-)S-} + 2{dx)xS+], 
jy = — [2cos^xWy+ + 2sm^ xuy- -sm2x{{uy+ + Uy^)C+ + {ux+ + Ux-)S-} + 2{dcc)xC-], 



with 



(3-1) 



«M± = Voi9t^'t>±-l\)- 



(3-2) 



In the stable and metastable states for given 9, the current perpendicular to the 
domain wall should vanish. If this were not the case, current would flow everywhere 
in the bulk, leading to a large energy cost. Neglecting the contributions of Ay leads 
to Eq. ( p-12 ) for A^, so that the variational solution already ensures that the current 
perpendicular to the wall vanishes for a of the (meta)stable state. 
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Fig. 3. Rectangular path encircling the flux line. We separate the path into parts 1, 2, 3 and 4 in 
order to calculation the contribution to the flux. 

We now consider vortices on the domain wall corresponding to the winding of 
the order parameter phase. By analogy to Josephson junctions, such a vortex may 
be considered a "kink" of the phase difference a. In the ordinary case, the size of 
the kink is an integer multiple of 27r. For the domain wall it is, however, possible to 
generate a kink to a metastable state. While, in general, a 27r-kink is associated with 
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the standard magnetic flux quantum ^o = hc/2e, a smaller kink between the stable 
and metastable state would carry a fraction of <Pq only, depending on the size of the 
kink and other parameters. For the Josephson junction we can describe the vortex 
using a Sine-Gordon equation of the phase difference containing only one length 
scale, the Josephson penetration depth. This analogy, however, is limited, because 
in the case of the domain wall two length scales have to be taken into account: the 
coherence length ^ (length scale of the variation of a) along the domain wall and the 
London penetration depth A^. We consider the limit Al S> ^, ignore the detailed 
structure of the vortex on small length scales (~ £,) and focus mainly on its magnetic 
properties. 

Let us analyze now a kink of a as a line defect in the domain wall between the 
stable state with a = uq and the metastable state with a = ai. We encircle this 
line by a wide rectangular path (see Fig. 3) and calculate the enclosed magnetic flux 
^ = § ds ■ A. First, we consider the parts of the path which cross the domain wall 
perpendicularly, 1 and 3. We choose the gauge such that dx(t)± = along both paths, 
i.e., any variation of (j)± is restricted to the parts of the path parallel to the domain 
wall. If we now use the condition jx = along this path, we obtain the contribution 



fi 



dxAx(x,v = —L) = / dB — 



^oniao,e) (3-3) 



to the flux from path 1, where 



n{a,9) 



S. 



27rJl - C2 



arctan 



C. 



'l-Cl 



Analogously, the contribution from path 3 is given by 



V's 



/ — CO 
dxAx{x,y = +L) = -<Pon{ai,0), 
,-oo 



(3-4) 



(3-5) 



and for the paths 2 and 4 we take the variation of (pj- into account. With jy = 
sufficiently far from the domain wall, this leads to 



-L 



f2 + f4 = dy{Ay{x ^ +oo,y) - Ay{x ^ -oo,y)) 
= - / (5y0+5j/0_) = ^0 



(3-6) 



Thus, the flux for the line defect depends on the angle 6 as 

ai — ao 



HG) = ^0 



27r 



+ iI(ao,6l)-77(ai, 



+ ^on, 



(3-7) 



which determines the flux up to an integer multiple of ^o- The possible magnetic 
fluxes are fractional and the smallest ones are smaller than $q in magnitude. This is 
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analogue Id tlie Jjpsephson vortices on a time-reversal symmetry breaking Josephson 
junction. E3*'I1ZP'I12P'Ej) In that case, however, the two junction states separated by the 
flux line are degenerate. This is the case here only for the special angles = ib7r/4. 
In the special case 9 = the phase difl'erences ckq = and ai = it lead to ^ = ±<Po/2, 
half a flux quantum. In all other cases the flux depends continuously on 6 and u. 

It is easy to see that two kinks in sequence (stable — > metastable — > stable) yield 
a total flux ^ = ^qu. Therefore, a conventional vortex (flux ^q) may decay into 
two fractional vortices on the domain wall (Fig. 4). Through this dissociation, the 
magnetic fleld energy can be reduced. If we assume that the domain wall is in the 
stable state, then the splitting of the vortex introduces a metastable domain wall 
connecting the two fractional flux lines. Thus, the cost of domain wall energy has to 
be compensated by the reduction of magnetic energy. 



^ 



<^ = <|)+<^ 



II <^ 



^ 



Fig. 4. A standard vortex with flux ^o on the domain wall decays into two fractional vortices whose 
fluxes add up to ^o- 

Let US compare the two energies here. The domain wall energy cost per unit 
area is edw = fdw{ai,0) - fdw{ao,0), with fdw defined in Eq. ( |2-17| ). The magnetic 
energy of the flux lines can be estimated by assuming that their fleld distribution is 
not much different from that of a standard vortex described by the London equation 
with the free energy functional 



Fl = J d^B^ + XliV X B)% 



(3-8) 



where we assume that the London penetration depth is not modifled by the domain 
wall. We now consider two fractional fluxes, ^1,^2 > 0, with ^1 + (^2 = ^0 (field 
along the z-axis) at the positions ri and r2, respectively. The resulting London 
equation is 



B, 



XlV^B, 



^i6^^\r-ri)+^2S'-^\r-r2), (3-9) 

which is taken in two dimensions, assuming homogeneity along the z-axis (r = (x, y)). 
The solution of the equation leads to the line energies ei and £2 for the two vortices 
and an interaction energy that depends on the distance between the vortices. Thus 
the total magnetic energy is 



v{\y2 -yi\) 






m -yi\ 



(3-10) 
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with Kq denoting the MacDonalds function. The hne energies contain as a lower 
cutoff length ^j (j = 1, 2), which is of a magnitude similar to the coherence length, 




. , sn^^M ~ I • (3-11) 

The potential for the fractional vortex at a distance R = \y2 — yi\ is therefore 

V{R)=v{R)+ed^R-eo, (3-12) 

where Eq is the line energy per unit length of a standard vortex. The energy loss 
due to the metastable state introduces a string potential between the two vortices. 
The optimal distance R is chosen to minimize the potential energy V{R), which 
corresponds to a pinning potential for the vortex. Note that this potential is only 
valid for R ^ S^j. Decaying into two fractional vortices leads to a very effective 
pinning, since the recombination is necessary in order to dislocate the vortex again 
from the domain wall. This occurs if the two fractional vortices are forced to approach 
to -R~|j. 

For a comparison of the energies involved, we estimate the magnetic energy 
gained by the decay into fractional vortices. 



<Pl 



(47rAL)2 



(3-13) 



where the London penetration depth is A| = l/87rKr/o7^. The maximal string energy 
is obtained for ^ = 0, where 

In this case the two energies become comparable for R just a fraction of Xl- For 
angles 9 closer to ±7r/4, the string potential is smaller and the optimal separation 
is larger. Furthermore, we find that increasing the anisotropy of the Fermi surface, 
denoted by i', tends to stabilize fractional vortices. The stability of the fractional 
vortices does also depend on temperature, because the internal structure of the 
domain wall can change with lowering temperature. B3) 

In this discussion we have assumed that the domain wall is an infinite plane with 
one fixed 9. In real materials this is not the case, because domain walls are pinned at 
impurities and lattices defects and, consequently, may change orientation by having 
"corners" . Since the phase value a for the stable domain wall structure depends on 
9, a kink of a should occur at every such corner and, according to our analysis, be 
accompanied by a flux line. Similarly, one can see that the crossing of two domain 
walls introduces a fractional vortex along the cutting line. In general, the phase and 
flux structure of a non-planar domain wall can be rather complicated. However, we 
do not consider these aspects any further here. 
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3.2. Barrier effect of a domain wall 

We now turn to the situation in which the superconductor is in the mixed phase 
and contains many vortices. If the conditions are appropriate, some of the vortices 
wiU be trapped by the domain wall and decay into fractional vortices. These flux 
lines are strongly pinned and may form a queue on the domain wall like the planks of 
a fence. Other vortices approaching the domain wall now are repelled by this vortex 
fence and cannot easily penetrate or traverse the domain wall (Fig. 5). In this way 
the domain wall acts as a very effective barrier for vortices. 

The question arises under which conditions the domain wall becomes permeable 
again. The vortices outside the domain wall generate a pressure on the fractional 
vortices, influencing their density and, hence, their mutual distance. With increasing 
external vortex density p = B/^q the distance R should shrink. This may be de- 
scribed qualitatively by adding a term Ocxt-R^ to the potential V{R) where acxt > 
grows with increasing p at the domain wall. Hence, it is clear that a growing den- 
sity of vortices close to the domain wall forces eventually the fractional vortices to 
recombine when i? ~ ^, the extension of the kinks of a. The recombined vortices 
are only weakly pinned so that the fence becomes permeable. This process needs 
not occur over all the domain wall simultaneously , but may start with local leaks 
for vortices to pass through. A simple picture of the "barrier height" of the domain 
wall can be obtained by the following argument. The density of fractional vortices 
is more or less equal to that of the vortices immediately outside of the domain wall. 



Therefore, the relation for the spacing i? ~ p 



-1/2 



{^o/B) ^''^ leads to the barrier 



field B* oc ^^^. The temperature dependence of ^ close to the transition temperature 
(the second transition in the case of a system with double transition) is essentially 
proportional to \T — Tc\~^''^ so that the barrier field is growing as the temperature 
is lowered B* oc \T — TJ. 






• • 



Fig. 5. Schematic picture of a domain wall occupied by fractional flux line (shaded) which repel 
approaching vortices (filled circles) so that they do not traverse the domain wall. 



With these properties, the presence of domain walls should have considerable 
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influence on the flux motion in the lower range of magnetic fields in the mixed phase. 
In particular this is true for phenomena like flux creep and hysteresis behavior of 
the mixed phase. We first discuss the influence of the barrier effect on the flux creep 
which can be observed in the relaxation of remnant magnetization after a magnetic 
field has been applied for some time to a superconductor and then switched off. 
The flux relaxation behavior has been investigated in much detail by the group 
of Mota fqrthe heavy Fermion superconductors UPts and Ui-ajThajBeis and for 
Sr2Ru04.B^'l5''Q^ The first two superconductors exhibit double transitions, where the 
low-temperature phase in each case is very likely time-reversal symmetry breaking 
and should therefore provide the conditions for domain walls as described above. 
Mota and coworkers observed that the initial flux relaxation rate drops drastically 
below the second superconducting transition. Hence, many vortices generating the 
remnant magnetization find it apparently more difficult to escape from the sample 
in the low-temperature phase. While initially almost no remnant flux leaves, the 
relaxation recovers after a longer waiting time. The absence of flux relaxation can 
be understood as the barrier effect of the domain wall. In the low-temperature 
phase, domains appear covering the sample with many domain walls. Vortices of 
the remnant magnetization are fenced in by these walls, and only a small portion 
of the total flux (mainly that close to surface not impeded by barriers) can move 
out when the external field is turned off. However, the encircled vortices press the 
domain walls which then may move slowly, and, finally, new pathways can open for 
the flux lines to leave the sample. Thus, the later appearance of flux decay can be 
attributed to slow domain wall motion. In the high-temperature phase domain walls 
are absent leading to flux creep as usual. 

The condition for this behavior is the existence of a large number of domain 
walls throughout the sample which are pinned at impurities and lattice defects such 
that they do not move too easily. Domains are believed to nucleate randomly at 
the second superconducting transition if there is no bias for one type of domain. If 
the superconducting state breaks time-reversal symmetry, an external magnetic field 
could provide a bias. Therefore we may expect that if the superconductor is cooled 
in a magnetic field (at least for the time-reversal symmetry breaking transition), 
the distribution of domains should be unbalanced in favor of one type reducing the 
number of domain walls. Then, the flux creep, previously impeded by domain walls, 
would be more ordinary. ElF A similar picture arises also from the measurement of 
flux motion in the low-temperature phase of Ui_2;Tha;Bei3 by Zieve et aL.EHP The 
samples which are field cooled and zero-field cooled exhibit a difference with regard 
to the motion of flux. In the former case, flux lines move more easily through the 
sample than in the latter, since in the zero-field cooled case, the presence of domain 
walls would also here impede the flux flow. Thus, there is a clear difference between 
the field- and zero-field cooled situation. The effect continuously grows below the 
second transition, indicating a gradual increase of the barrier strength as suggested 
above. This is ari aspect that also appears in the hysteresis experiment discussed in 
the next section. B 
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§4. Intrinsic noise in the hysteresis 

The hysteresis of the magnetization in the mixed phase provides one way to 
measure the critical current. The basic features of the hysteresis cycle are described 
well by the critical voxtex state models among which the Bean model is the most 
simple representation. E3 In this model the vortices are described as a magnetic flux 
density. We omit here a detailed introduction, since the basic properties of the Bean 
model can be found in many textbooks. eP 




(a) 
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Fig. 6. Standard Bean profile of tlie flux density for the magnetization process in the hysteresis of 
a superconductor in the mixed state. The magnetization M(H) is the integral of flux density 
in the sample. 

For the following analysis we use the model of a superconducting slab which 
has a finite width along the x-axis and is infinite along the other two axes. Let 
us first study the standard behavior of the hysteresis, including the surface barrier 
effect of the lower critical field (or alternatively the Bean-Livingston barrier). We 
start by assuming that the external field was first raised to the maximum value Hm- 




Fig. 7. Standard hysteresis in the mixed state of a superconductor based on the Bean critical state 
model. 
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Then the field distribution is similar to that in Fig. 6a). Now we turn the field 
oflF gradually and reverse it until we reach —Hm (Fig. 6f). On the way the field 
distribution passes through the distributions shown in Fig. 6b)-e). It is important 
that in order to introduce reversed vortices into the system, the external field must 
exceed the lower critical field. Hence, in the range between -ffext = and —Hd, 
the internal flux distribution is unchanged, and at —Hd the magnetization jumps 
abruptly. After reaching i/ext = —Hm, the field is again gradually turned to zero and 
further to Hm, so that the hysteresis cycle closes. The magnetization M as a function 
of H corresponds to the integrated flux density in the Bean critical state. The slope 
of the fleld distribution corresponds to the critical current {Jc = {c/ATr)dxBz) and is 
proportional to the width of the hysteresis loop (Fig. 7) . 

If we now introduce domain walls as additional barriers within the sample, new 
structures appear in the hysteresis cycle. The domain wall is characterized by a 
speciflc barrier field (height) B*; i.e., if the local field at the domain wall exceeds 
B*, the domain wall is permeable to vortices. For simplicity we consider the case of 
two domain walls which fence in some region inside the slab. We start again from 
the situation in which the external field is at the value Hm, so that everywhere in 
the sample the local field value is larger than B* , and the Bean profile has developed 
everywhere with a slope corresponding to the critical current due to usual pinning 
(Fig. 8a). Now, we gradually reduce the external field. As we reach zero, we find 
that the domain walls have trapped some fiux; that is, more flux than usual is 
remaining in the sample (Fig. 8b). If we now turn to negative flelds, then initially 
at —Hd, reversed flux lines enter and annihilate the flux close to surface (Fig. 8c). 
With decreasing external fleld, the reversed flux lines reach the domain walls and 
annihilate vortices on the domain wall, which are in turn immediately replaced by 
new positive flux lines from the interior of the fenced in region. This process continues 
until the density of positive vortices on the domain wall is essentially zero and leads 
to a sharp increase in the magnetization similar to that at the lower critical field 
(Fig. 8d). Note that the negative vortices from the outside cannot enter the inner 
region yet, since at the end of this process the domain wall is occupied by reversed 
(negative) fractional vortices, which now provide a barrier. The fiux density must be 
enhanced from outside, and the domain wall only becomes permeable when the local 
field of the negative vortices has reached —B* (Fig. 8e). Then, negative vortices can 
pass freely through the domain wall and annihilate the remaining positive vortices 
on the other side. This leads to a further avalanche-like rise in the magnetization. 
Then the external field is lowered until we reach —Hm (Fig. 8f). The analogous 
continuation for increasing field finally leads to a closed hysteresis cycle. 

The additional sharp avalanche-like structures in the magnetization depend 
strongly on the position of the domain walls. As mentioned earlier, the domain 
walls are not rigidly fixed, but only pinned, and can change their position. During 
several cycles, the domain walls may be located at different places from cycle to 
cycle due to pressure from the vortices. Hence, the hysteresis cycles do not retrace 
themselves completely, but rather have strong deviations in the field range where 
the avalanche effects occur, i.e., close to the lower critical field, ±Hd- Taking the 
average over several cycles, the standard deviation, 6M = {{M{H) — {M{H)))'^) re- 



16 



M. Sigrist and D. F. Agterberg 




(b) 



(c) 



- 

-B * 



(d) 



(e) 




B 

-Hcl 



Fig. 8. Bean profile of a superconducting slab with two domain walls in a magnetic field, 
critical field and the barrier field are indicated on the axis on the right hand side. 
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Fig. 9. Hysteresis of a superconducting slab containing domain walls as barriers. Upper figure: 
Several cycles with varying location of the domain wall. Lower figure; Average of several cycles 
and standard deviation, which possess maxima around the region of ±Hci- 



veals a strong "noise" signal, as shown in Fig. 9 ((...) denotes the average over many 
cycles). This noise is intrinsic to the superconducting state with internal barriers. It 
is easy to see that the magnitude of this noise depends on the barrier height given by 
B* , and it should disappear, if B* goes to zero (the case of a completely permeable 
domain wall): 

5M oc B*5x, (4-1) 

where 5x is the standard deviation of the domain wall positions, only weakly depen- 
dent on temperature. 

We may compare this behavior now to the experiments done on UPts by Shung 
et aL.U^ The features of the noise in the hysteresis are qualitatively very similar to 
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those obtained in our simple simulation. In the high-temperature phase, the noise 
in the hysteresis cycle is small and featureless. However, simultaneously with the 
onset of the second transition, the intrinsic noise appears continuously and grows 



essentially linearly with \T — Tc2\ as we expect from Eq. (4-1). Within the cycle 
the noise appears in the same field range and with a similar structure, as in our 
simulation, i.e. around Hd- This suggests that the continuous increase of the noise 
is very likely connected with the broken time-reversal symmetry below the second 
transition where the degeneracy of the superconducting phase allows for domain 
formation. The barrier height introduced by these domain walls increases as the 
order parameter of the low-temperature phase grows. A similar picture arises also 
from the measurement of fiux motion in the low-temperature phase of Ui-^iTh^^Beis 
by Zieve et al.,B3) as mentioned above. 

§5. Conclusion 

We have shown that domain walls in time-reversal symmetry-breaking super- 
conductors can play an important role for the motion of flux lines. Under certain 
conditions, domain walls can accommodate flux lines whose flux is a fraction of the 
standard flux quantum. On such domain walls, conventional vortices can decay into 
two fractional vortices. If such fractional vortices line up, the domain wall becomes 
a very efficient barrier to vortices. 

Evidence for this property of the domain walljs eiven bv the absence of flux 
creep or in the intrinsic noise of the hysteresis cycle. IHP'u)'q)'H)'E3) In the case of UPta 
and Ui_2^Th2.Bei3 experiments reveal a clear connection between the onset of the 
time-reversal symmetry-breaking in the low-temperature phase and the appearance 
of a new flux trapping mechanism. In Sr2Ru04 the drop of the creep rate is not 
connected with the onset of the time-reversal symmetry breaking phase. Rather the 
drop around 50 mK may be associated with a "transition" of the domain wall states. 
This may be associated with the multi-band nature of the superconducting state, 
which can lead to the temperature dependence of parameters like the anisotropy 
factor V, since the relative contribution to superconductivity by the different bands 
certainly varies with temperature. This case needs definitely further consideration. 
The barrier effect of domain walls does not, in general, influence the magnitude 
of the critical current Jc, which is related to ordinary impurity- induced pinning 
effects. Although the initial flux relaxation is vanishing in the low-temperature 
phase, Jc remains flnite and does not show any anomalous temperature dependence, 
indicating that the phenomena are not connected with a change in the pinning of 
individual vortices. We may take the observation of such flux flow phenomena as 
indirect evidence for the presence of domain walls. To this time in none of the 
above mentioned superconductors has the direct observation of the domain walls 
been reported. Clearly, the discovery of lined-up fractional vortices could be one 
direct experimental veriflcation. Another proposal is based on the modification of 
the local quasiparticle density of states in the domain wall which may be observable 
by scanning tunneling microscopy. EJ) This paper presents a plausible explanation for 
the basic properties seen in the experiments of three unconventional superconductors. 
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Nevertheless, several details in the experimental data, which may be typical for the 
phenomenon or specific to each system, remain to be discussed in future. 
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